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Wo have searched for a heavy resonance decaying into a Z+]ct final state in pp collisions at a center 
of mass energy of 1.96 TeV at the Fermilab Tevatron collider using the DO detector. No indication for 
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such a resonance was found in a data sample corresponding to an integrated luminosity of 370 pb~ . 
We set upper limits on the cross section times branching fraction for heavy resonance production at 
the 95% C.L. as a function of the resonance mass and width. The limits are interpreted within the 
framework of a specific model of excited quark production. 

PACS numbers: 13.85.Rm,14.65.-q,14.70.Hp,14.80.-j 



Heavy resonances decaying into a quark and a gauge 
boson may signal the existence of excited quarks and 
thereby indicate quark substructure 51 ■ Searches for ex- 
cited quarks have been carried out in the past using dijet 
[llEHi photon+jet, and W+jet final states. In the 
analysis described here, we searched for resonances in 
the Z+jet channel, where the Z boson is detected via its 
Z e+e^ decay mode. This signature is practically free 
of instrumental background. However, it suffers from the 
low branching fraction (3.36%) of the Z — > e+e~ decay 
channel. The high luminosity delivered by the Fermilab 
Tevatron collider in Run H makes it possible to present 
results on this final state for the first time. 

For the production and decay of a resonance, we con- 
sidered the model [li implemented in PYTHIA 6.202 |^. 
Here, a quark (antiquark) and a gluon from the colliding 
proton and antiproton form a resonance, q*, which subse- 
quently decays into a Z boson and a quark: — > g -I- Z. 
The model has two free parameters, Mq* , the mass of 
the resonance, and A, the compositeness scale. They 
determine the production cross section and the natural 
width of the resonance. The latter scales as 1/^^, where 

The Run II DO detector Q consists of several layered 
subdetectors. For the present analysis, the most rele- 
vant parts are the liquid- argon/uranium calorimeter Q 
and the central tracking system. The calorimeter, di- 
vided into electromagnetic and hadronic sections, has a 
granularity of A77 x = 0.1 x 0.1, where t] is the pseu- 
dorapidity (r\ = — In [tan(0/2)] with 6 being the polar 
angle measured from the geometrical center of the de- 
tector with respect to the proton beam direction) and 
(j) is the azimuthal angle. The third innermost layer, in 
which the largest electromagnetic energy deposition is ex- 
pected, has a finer granularity of Ar/ x A^ — 0.05 x 0.05. 
The central calorimeter covers \r]\ < 1.1, and the two end 
calorimeters extend coverage to \ri\ « 4.5. The tracking 
system consists of a silicon microstrip tracker and a cen- 
tral fiber tracker, both located within a 2 T supercon- 
ducting solenoidal magnet, with designs optimized for 
tracking and vertexing at pseudorapidities I77I < 3 and 
1 77 1 < 2, respectively. 

The data used in this analysis were collected between 
April 2002 and August 2004, with an integrated lumi- 
nosity of 370 pb^^. The selected events were required to 
pass at least one of several single- or di-electron triggers. 
The efficiency of the trigger was measured with data and 
found to reach a plateau of etrig — 0.982±0.011 for events 
satisfying the final event selection criteria. 



Offline event selection was based on run quality, event 
properties, and electron and jet identification criteria. 
Events were required to have a reconstructed vertex with 
a longitudinal position within 60 cm of the detector cen- 
ter. Electrons were reconstructed from electromagnetic 
(EM) clusters in the calorimeter using a cone algorithm. 
The reconstructed electron candidates were required to 
satisfy either \ri\ < 1.1 or 1.5 < I77I < 2.5. Electron pairs 
with > 30 GeV and p§? > 25 GeV in the event were 
used to reconstruct the Z boson candidate. The electron 
pair was required to have an invariant mass M^e near the 
Z boson mass, 80 < Mee < 102 GeV. 

To reduce background contamination, mainly from jets 
misidentified as electrons, the EM clusters were required 
to pass three quality criteria based on shower profile: (z) 
the ratio of the energy deposited in the electromagnetic 
part of the calorimeter to the total shower energy had 
to exceed 0.9; (ii) the lateral and longitudinal shapes 
of the energy cluster had to be consistent with those 
of an electron; and (iii) the electron had to be isolated 
from other energy deposits in the calorimeter with iso- 
lation fraction fiso < 0.15. The isolation fraction is de- 
fined as fiso = [E{OA) - Eem{0-2)] /Eem{0.2), where 
E{Rconc) and Eem (Rconc) are the total and the EM 
energy, respectively, deposited within a cone of radius 
-Rconc — \J (A?7)2 -I- (A(/))2 centered around the electron. 
Additionally, at least one of the electrons was required 
to have a spatially close track with a momentum con- 
sistent with the EM shower energy. A total of 24,734 
events passed these criteria. In Fig. ^ the distribution 
of the invariant mass, Mee, of the two selected electrons 
is shown. A very clean, almost background-free Z boson 
signal is evident. 

Jets were reconstructed using the "Run II cone algo- 
rithm" 9] which combines cell energies within a cone of 
radius i?conc — 0.5. Spurious jets from isolated noisy 
calorimeter cells were supressed by cuts on the jet shape 
and by requiring that the charged tracks associated with 
the jet had to carry a minimum fraction of the jet trans- 
verse energy. The transverse momentum of each jet was 
corrected for offsets due to the underlying event, multiple 
interactions and noise, out-of-cone showering, and the 
detector energy response as determined from the trans- 
verse energy balance of photon-|-jet events. Jets were 
required to have pt > 20 GeV and jryj < 2.5 and to not 
overlap with any of the reconstructed EM objects within 
a distance of 0.4 in (77, 0) space. Requiring one or more 
jets with these selection criteria, 2,417 data events re- 
main. 
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FIG. 1: The invariant mass of the two selected electrons in 
the data events. 

We have considered two kinds of instrumental back- 
grounds where hadronic jets are misreconstructed as EM 
clusters and mimic Z boson events. A background from 
genuine QCD multi-jet production arises when both of 
the EM objects are hadronic jets that fluctuate to elec- 
tromagnetic final states. This background has been es- 
timated to be (0.56 ± 0.02)% of the signal in the mass 
region of 80 < Al^e < 102 GeV as calculated by com- 
paring the Mee mass distribution of the selected events 
with a distribution that required inverted shower shape 
criteria and the absence of matching tracks. The other 
source of background is Vl^ — s- ei^-|-jets events where a 
hadronic jet is misidentified as an electron. These events 
are characterized by significant missing transverse energy 
(^t), and should also appear in the data sample where 
only one of the EM objects has a matched track. From 
comparison of the distribution of these events with 
that where both electrons do have matched tracks, we 
estimate that this background is an order of magnitude 
less than the QCD background. 

The main standard model (SM) background to the ex- 
cited quark signal is inclusive Z/j* — s- e+e^ pair pro- 
duction which has been simulated with pythia using the 
CTEQ5L 113 parton distribution functions (PDFs). In 
order to enhance the statistics for events where the in- 
variant mass of the Z boson and the leading jet, Mzji, 
is high, in addition to the so-called 2^1 process, we 
have also generated events including matrix elements of 
first order in (2 — > 2 process) with different thresholds 
of Mzq, the invariant mass of the Z boson and the ac- 
companying parton in the final state. A minimum value 
of 30 GeV for pTp, the transverse momentum of the par- 
ton in the 2-^2 collision, has been set in order to avoid 



FIG. 2: Invariant mass distribution of the Z boson and the 
leading jet, Mzji- The data are shown by the full squares 
with error bars. The actual number of events in a bin is the 
product of the plotted value and the bin width measured in 
10 GeV units. The SM backgrounds generated with pythia 
are shown in the histograms: 2 — > 1 without threshold (solid 
line) ,2^2 with various Mzq thresholds (discontinuous lines, 
as indicated). Each curve with a defininte Alzq threshold 
value stops when the curve of the next threshold value takes 
over. Also shown with open circles is the signal due to an 
excited quark of 500 GeV mass and narrow width (^ = 1). 
The resonance production cross section is taken from Ref. 

Q. 



collinear divergences. The leading jet pt distribution has 
a mean and an RMS value of 106 GeV and 27 GeV, re- 
spectively, at the lowest resonance mass investigated and 
after the final selection, therefore the pTp cut does not 
affect the analysis. The shape of the Mzji distribution 
has been compared with that obtained with the ALPGEN 
program 

'ill 

and there is reasonable agreement between 
them. Any differences in the background level have been 
taken into account as a systematic uncertainty. 

Signal events were generated with PYTHIA using the 
CTEQ5L PDFs for the following resonance mass val- 
ues: Mg. = 300, 400, 500, 600 and 700 GeV with ^ 
= K/Mq- = 1. For each mass, except for the lowest 
one, we also generated events with ^ = 0.3, 0.5, 0.7, in 
order to vary the natural width of the resonance, Vq- . 
The form factors associated with the interaction of the 
quarks with the SM gauge bosons were set to unity. The 
MC events were passed through the same reconstruction 
software and selection criteria as the data. The events 
have been used to estimate the geometrical acceptance 
and jet and electron identification efficiencies. The com- 
bined acceptance times efficiencies are listed in Table 
The resolution of Mzji has been found to be ~ 9%. 
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In Fig.[21we compare the Mzji distribution of the data 
with the PYTHIA 2^1 process and with the pythia 
2^2 processes with various Mzq thresholds. For the 
2^1 process, the MC is normahzed to the total number 
of data events. A different but common normalization 
factor is used for all 2 ^ 2 processes determined using 
the Mzq > 100 GeV MC sample for Mzji > 150 GeV. 
The 2^1 simulation agrees well with the data but pro- 
vides sufficient statistics only for Mz ji < 300 GeV. On 
the other hand, the 2 — > 2 processes describe the data 
with reasonable precision for Mzji > 150 GeV. Since 
the latter is the region of interest for the present search, 
we have used only the 2 — > 2 process for estimation of the 
SM background with an Mzq threshold chosen according 
to the Mzji region to be investigated. Also shown in 
Fig. 13 is the signal due to an excited quark of 500 GeV 
mass and narrow width — I). 



400: 
350: 
300: 
250: 

> 

(u : 
2,200: 



Mq.=500 GeV 
SM background 




„. .-■•■^.VA^s.'-;:-:;. 




°0 100 200 300 400 500 600 700 800 
Mzji [GeV] 

FIG. 3: pTZ vs Mzji distributions for a resonance of mass of 
500 GeV = 1) and for the SM background. Both the signal 
and background events passed through complete reconstruc- 
tion. Each distribution is arbitrarily normalized. 



Since no significant excess of events is observed, which 
would indicate the presence of a resonance, we deter- 
mined the upper limit on the production cross section of 
a hypothetical resonance as a function of its mass and 
width. We made use of the fact that in the prz vs Mzji 
plane, events from the resonance are concentrated for 
Mzji around the mass value and for prz at about half 
of the mass value of the resonance, since the resonance 
is nearly at rest. The SM background does not exhibit a 
similar structure, as it is shown in Fig. 13 In addition, fi- 
nite width and mass resolutions wash out the correlation 
between prz and Mzji- We therefore considered events 
around the peak values M^j^ and of the resonance 



determined by the following condition: 
f Mzji 



Mrms 



PTZ - PtZ 

rms 

Ptz 



(1) 



and we optimized the cut value k. Here, and p^^'^ 

are the RMS values of the corresponding distributions of 
the resonance. At given values of mass and width, the 
latter defined by ^, we varied k in Eqn.l^ between and 
3 in steps of 0.1. 



TABLE I: Measured (ergs) and expected (0-95*) values of the 
upper limit on the resonance cross section times branching 
fraction, signal acceptance x efficiency, SM background, and 
number of observed events at the optimal value of the topo- 
logical cut k for different resonance masses and for ^ = 1. 





k 


0"95 


ave 


Acceptance 


SM 


Data 


(GeV) 




(Pb) 


(pb) 


X efficiency 


background (events) 


300 


1.1 


0.25 


0.290 


0.140 ± 0.009 


32.8 ± 2.9 


31 


400 


1.2 


0.15 


0.129 


0.164 ± 0.010 


7.5 ± 0.8 


9 


500 


1.3 


0.08 


0.079 


0.195 ± 0.012 


2.9 ± 0.8 


3 


600 


1.8 


0.05 


0.053 


0.244 ± 0.014 


1.6 ± 0.6 


1 


700 


1.7 


0.03 


0.044 


0.243 ± 0.014 


0.64 ± 0.06 






TABLE II: Measured upper limit on the resonance cross sec- 
tion times branching fraction at the 95% C.L., (T95, for differ- 
ent resonance masses and ^ values, aq* , the production cross 
section of an excited quark times its decay branching fraction 
into Z+jet and Z — > e^e~, is calculated in LO. The width 
for ^ = 1 is also shown Q for each resonance mass. Cross 
sections are quoted in pb, whereas masses and widths are in 
GeV. 



Mass 




0"95 






0-5. 






C = 0.3 


0.5 


0.7 


1 




1 


300 








0.25 


2.045 


13 


400 


0.32 


0.16 


0.15 


0.15 


0.382 


16 


500 


0.17 


0.08 


0.07 


0.08 


0.084 


20 


600 


0.10 


0.06 


0.05 


0.05 


0.021 


24 


700 


0.07 


0.05 


0.05 


0.03 


0.005 


27 



Based only on the information from the signal and 
background simulation, for each k we calculated crgg*^, 
the expected value of the upper limit on the resonance 
production cross section times branching fraction at the 
95% C.L. using a Bayesian approach [13 and by aver- 
aging over possible outcomes of the background-only hy- 
pothesis assuming Poisson statistics of the background. 
The optimum value of k corresponds to the minimum 



value of crgg*^. At this value of fc, using also the data, we 
derived 1795 , the measured value of the upper limit on the 
resonance production cross section times branching frac- 
tion at the 95% C.L. In this calculation we have taken 
into account systematic uncertainties in the determina- 
tion of the luminosity (6.5%), trigger and identification 
efficiencies, and those of the jet calibration and resolu- 
tion. Systematic uncertainties due to the modeling of 
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FIG. 4: Upper limit on the resonance cross section times 
branching fraction at the 95% C.L., ugs, for different reso- 
nance masses as a function of the resonance width. 
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FIG. 5: Upper limits on the resonance cross section times 
branching fraction at the 95% C.L., ergs, for different ^ values 
as functions of the resonance mass (open symbols). Full circles 
indicate the LO production cross section of an excited quark 
times its decay branching fraction into Z+}et and Z — > e'^e~ , 
ag„ for C = 1 Q- 



the SM background and to the choice of the PDF, as 
well as those due to the threshold of the Mzji > 150 
GeV in the normalization of the background, have also 
been included. 

In Table m 0-95 and ag^'^ are shown together with the 
signal acceptance, the SM background level, and the 



number of data events for ^ = 1. The measured ergs 
values are displayed in Fig. ^ and Fig. [SJ and are com- 
piled in Table Ull for different masses and widths {£,). In 
Fig. 13 also shown is cr^* , the LO production cross sec- 
tion of an excited quark times its decay branching frac- 
tion into Z-l-jet and Z e+e^, for <^ = 1 ,lj. We find 
a lower limit of 510 GeV at the 95% C.L. for the mass 
of an excited quark for ^ = 1 within the framework of 
the model considered. In earlier measurements, lower 
bounds of 460, 530 and 775 GeV were obtained for the 
same quantity, but in different decay modes, namely in 
q* 57 5], q* qW 0, and 9* ^ (?5 0], respectively, 
and therefore with different systematics. 

In conclusion, we have searched for a resonance pro- 
duced by the fusion of a gluon and a quark in pp collisions 
at a center of mass energy of 1.96 TeV which decays into 
a Z boson and a quark in the Z —^ e^e~ decay chan- 
nel. In the absence of a signal, we have determined 95% 
C.L. upper limits on the cross section times branching 
fraction as a function of the mass and width of the res- 
onance. The present study is complementary to earlier 
searches because it has sensitivity to hypothetical models 
with enhanced couplings to the Z boson. 
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